\(\int \sec (c+d x) (a+a \sec (c+d x))^{5/2} (A+C \sec ^2(c+d x)) \, dx\) [175]
Optimal result
Integrand size = 33, antiderivative size = 169 \[
\int \sec (c+d x) (a+a \sec (c+d x))^{5/2} \left (A+C \sec ^2(c+d x)\right ) \, dx=\frac {64 a^3 (21 A+13 C) \tan (c+d x)}{315 d \sqrt {a+a \sec (c+d x)}}+\frac {16 a^2 (21 A+13 C) \sqrt {a+a \sec (c+d x)} \tan (c+d x)}{315 d}+\frac {2 a (21 A+13 C) (a+a \sec (c+d x))^{3/2} \tan (c+d x)}{105 d}-\frac {4 C (a+a \sec (c+d x))^{5/2} \tan (c+d x)}{63 d}+\frac {2 C (a+a \sec (c+d x))^{7/2} \tan (c+d x)}{9 a d}
\]
[Out]
2/105*a*(21*A+13*C)*(a+a*sec(d*x+c))^(3/2)*tan(d*x+c)/d-4/63*C*(a+a*sec(d*x+c))^(5/2)*tan(d*x+c)/d+2/9*C*(a+a*
sec(d*x+c))^(7/2)*tan(d*x+c)/a/d+64/315*a^3*(21*A+13*C)*tan(d*x+c)/d/(a+a*sec(d*x+c))^(1/2)+16/315*a^2*(21*A+1
3*C)*(a+a*sec(d*x+c))^(1/2)*tan(d*x+c)/d
Rubi [A] (verified)
Time = 0.40 (sec) , antiderivative size = 169, normalized size of antiderivative = 1.00, number of
steps used = 5, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.121, Rules used = {4168, 4086, 3878, 3877}
\[
\int \sec (c+d x) (a+a \sec (c+d x))^{5/2} \left (A+C \sec ^2(c+d x)\right ) \, dx=\frac {64 a^3 (21 A+13 C) \tan (c+d x)}{315 d \sqrt {a \sec (c+d x)+a}}+\frac {16 a^2 (21 A+13 C) \tan (c+d x) \sqrt {a \sec (c+d x)+a}}{315 d}+\frac {2 a (21 A+13 C) \tan (c+d x) (a \sec (c+d x)+a)^{3/2}}{105 d}+\frac {2 C \tan (c+d x) (a \sec (c+d x)+a)^{7/2}}{9 a d}-\frac {4 C \tan (c+d x) (a \sec (c+d x)+a)^{5/2}}{63 d}
\]
[In]
Int[Sec[c + d*x]*(a + a*Sec[c + d*x])^(5/2)*(A + C*Sec[c + d*x]^2),x]
[Out]
(64*a^3*(21*A + 13*C)*Tan[c + d*x])/(315*d*Sqrt[a + a*Sec[c + d*x]]) + (16*a^2*(21*A + 13*C)*Sqrt[a + a*Sec[c
+ d*x]]*Tan[c + d*x])/(315*d) + (2*a*(21*A + 13*C)*(a + a*Sec[c + d*x])^(3/2)*Tan[c + d*x])/(105*d) - (4*C*(a
+ a*Sec[c + d*x])^(5/2)*Tan[c + d*x])/(63*d) + (2*C*(a + a*Sec[c + d*x])^(7/2)*Tan[c + d*x])/(9*a*d)
Rule 3877
Int[csc[(e_.) + (f_.)*(x_)]*Sqrt[csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_)], x_Symbol] :> Simp[-2*b*(Cot[e + f*x]/(
f*Sqrt[a + b*Csc[e + f*x]])), x] /; FreeQ[{a, b, e, f}, x] && EqQ[a^2 - b^2, 0]
Rule 3878
Int[csc[(e_.) + (f_.)*(x_)]*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_), x_Symbol] :> Simp[(-b)*Cot[e + f*x]*(
(a + b*Csc[e + f*x])^(m - 1)/(f*m)), x] + Dist[a*((2*m - 1)/m), Int[Csc[e + f*x]*(a + b*Csc[e + f*x])^(m - 1),
x], x] /; FreeQ[{a, b, e, f}, x] && EqQ[a^2 - b^2, 0] && GtQ[m, 1/2] && IntegerQ[2*m]
Rule 4086
Int[csc[(e_.) + (f_.)*(x_)]*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_)*(csc[(e_.) + (f_.)*(x_)]*(B_.) + (A_))
, x_Symbol] :> Simp[(-B)*Cot[e + f*x]*((a + b*Csc[e + f*x])^m/(f*(m + 1))), x] + Dist[(a*B*m + A*b*(m + 1))/(b
*(m + 1)), Int[Csc[e + f*x]*(a + b*Csc[e + f*x])^m, x], x] /; FreeQ[{a, b, A, B, e, f, m}, x] && NeQ[A*b - a*B
, 0] && EqQ[a^2 - b^2, 0] && NeQ[a*B*m + A*b*(m + 1), 0] && !LtQ[m, -2^(-1)]
Rule 4168
Int[csc[(e_.) + (f_.)*(x_)]*((A_.) + csc[(e_.) + (f_.)*(x_)]^2*(C_.))*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(
m_), x_Symbol] :> Simp[(-C)*Cot[e + f*x]*((a + b*Csc[e + f*x])^(m + 1)/(b*f*(m + 2))), x] + Dist[1/(b*(m + 2))
, Int[Csc[e + f*x]*(a + b*Csc[e + f*x])^m*Simp[b*A*(m + 2) + b*C*(m + 1) - a*C*Csc[e + f*x], x], x], x] /; Fre
eQ[{a, b, e, f, A, C, m}, x] && !LtQ[m, -1]
Rubi steps \begin{align*}
\text {integral}& = \frac {2 C (a+a \sec (c+d x))^{7/2} \tan (c+d x)}{9 a d}+\frac {2 \int \sec (c+d x) (a+a \sec (c+d x))^{5/2} \left (\frac {1}{2} a (9 A+7 C)-a C \sec (c+d x)\right ) \, dx}{9 a} \\ & = -\frac {4 C (a+a \sec (c+d x))^{5/2} \tan (c+d x)}{63 d}+\frac {2 C (a+a \sec (c+d x))^{7/2} \tan (c+d x)}{9 a d}+\frac {1}{21} (21 A+13 C) \int \sec (c+d x) (a+a \sec (c+d x))^{5/2} \, dx \\ & = \frac {2 a (21 A+13 C) (a+a \sec (c+d x))^{3/2} \tan (c+d x)}{105 d}-\frac {4 C (a+a \sec (c+d x))^{5/2} \tan (c+d x)}{63 d}+\frac {2 C (a+a \sec (c+d x))^{7/2} \tan (c+d x)}{9 a d}+\frac {1}{105} (8 a (21 A+13 C)) \int \sec (c+d x) (a+a \sec (c+d x))^{3/2} \, dx \\ & = \frac {16 a^2 (21 A+13 C) \sqrt {a+a \sec (c+d x)} \tan (c+d x)}{315 d}+\frac {2 a (21 A+13 C) (a+a \sec (c+d x))^{3/2} \tan (c+d x)}{105 d}-\frac {4 C (a+a \sec (c+d x))^{5/2} \tan (c+d x)}{63 d}+\frac {2 C (a+a \sec (c+d x))^{7/2} \tan (c+d x)}{9 a d}+\frac {1}{315} \left (32 a^2 (21 A+13 C)\right ) \int \sec (c+d x) \sqrt {a+a \sec (c+d x)} \, dx \\ & = \frac {64 a^3 (21 A+13 C) \tan (c+d x)}{315 d \sqrt {a+a \sec (c+d x)}}+\frac {16 a^2 (21 A+13 C) \sqrt {a+a \sec (c+d x)} \tan (c+d x)}{315 d}+\frac {2 a (21 A+13 C) (a+a \sec (c+d x))^{3/2} \tan (c+d x)}{105 d}-\frac {4 C (a+a \sec (c+d x))^{5/2} \tan (c+d x)}{63 d}+\frac {2 C (a+a \sec (c+d x))^{7/2} \tan (c+d x)}{9 a d} \\
\end{align*}
Mathematica [A] (verified)
Time = 2.26 (sec) , antiderivative size = 90, normalized size of antiderivative = 0.53
\[
\int \sec (c+d x) (a+a \sec (c+d x))^{5/2} \left (A+C \sec ^2(c+d x)\right ) \, dx=\frac {2 a^3 \left (903 A+584 C+(294 A+292 C) \sec (c+d x)+3 (21 A+73 C) \sec ^2(c+d x)+130 C \sec ^3(c+d x)+35 C \sec ^4(c+d x)\right ) \tan (c+d x)}{315 d \sqrt {a (1+\sec (c+d x))}}
\]
[In]
Integrate[Sec[c + d*x]*(a + a*Sec[c + d*x])^(5/2)*(A + C*Sec[c + d*x]^2),x]
[Out]
(2*a^3*(903*A + 584*C + (294*A + 292*C)*Sec[c + d*x] + 3*(21*A + 73*C)*Sec[c + d*x]^2 + 130*C*Sec[c + d*x]^3 +
35*C*Sec[c + d*x]^4)*Tan[c + d*x])/(315*d*Sqrt[a*(1 + Sec[c + d*x])])
Maple [A] (verified)
Time = 44.86 (sec) , antiderivative size = 124, normalized size of antiderivative =
0.73
| | |
| method | result | size |
| | |
| default |
\(\frac {2 a^{2} \left (903 A \cos \left (d x +c \right )^{4}+584 C \cos \left (d x +c \right )^{4}+294 A \cos \left (d x +c \right )^{3}+292 C \cos \left (d x +c \right )^{3}+63 A \cos \left (d x +c \right )^{2}+219 C \cos \left (d x +c \right )^{2}+130 C \cos \left (d x +c \right )+35 C \right ) \sqrt {a \left (1+\sec \left (d x +c \right )\right )}\, \tan \left (d x +c \right ) \sec \left (d x +c \right )^{3}}{315 d \left (\cos \left (d x +c \right )+1\right )}\) |
\(124\) |
| parts |
\(\frac {2 A \,a^{2} \sqrt {a \left (1+\sec \left (d x +c \right )\right )}\, \left (43 \sin \left (d x +c \right )+14 \tan \left (d x +c \right )+3 \sec \left (d x +c \right ) \tan \left (d x +c \right )\right )}{15 d \left (\cos \left (d x +c \right )+1\right )}+\frac {2 C \,a^{2} \left (584 \cos \left (d x +c \right )^{4}+292 \cos \left (d x +c \right )^{3}+219 \cos \left (d x +c \right )^{2}+130 \cos \left (d x +c \right )+35\right ) \sqrt {a \left (1+\sec \left (d x +c \right )\right )}\, \tan \left (d x +c \right ) \sec \left (d x +c \right )^{3}}{315 d \left (\cos \left (d x +c \right )+1\right )}\) |
\(149\) |
| | |
|
|
|
|
[In]
int(sec(d*x+c)*(a+a*sec(d*x+c))^(5/2)*(A+C*sec(d*x+c)^2),x,method=_RETURNVERBOSE)
[Out]
2/315*a^2/d*(903*A*cos(d*x+c)^4+584*C*cos(d*x+c)^4+294*A*cos(d*x+c)^3+292*C*cos(d*x+c)^3+63*A*cos(d*x+c)^2+219
*C*cos(d*x+c)^2+130*C*cos(d*x+c)+35*C)*(a*(1+sec(d*x+c)))^(1/2)/(cos(d*x+c)+1)*tan(d*x+c)*sec(d*x+c)^3
Fricas [A] (verification not implemented)
none
Time = 0.26 (sec) , antiderivative size = 130, normalized size of antiderivative = 0.77
\[
\int \sec (c+d x) (a+a \sec (c+d x))^{5/2} \left (A+C \sec ^2(c+d x)\right ) \, dx=\frac {2 \, {\left ({\left (903 \, A + 584 \, C\right )} a^{2} \cos \left (d x + c\right )^{4} + 2 \, {\left (147 \, A + 146 \, C\right )} a^{2} \cos \left (d x + c\right )^{3} + 3 \, {\left (21 \, A + 73 \, C\right )} a^{2} \cos \left (d x + c\right )^{2} + 130 \, C a^{2} \cos \left (d x + c\right ) + 35 \, C a^{2}\right )} \sqrt {\frac {a \cos \left (d x + c\right ) + a}{\cos \left (d x + c\right )}} \sin \left (d x + c\right )}{315 \, {\left (d \cos \left (d x + c\right )^{5} + d \cos \left (d x + c\right )^{4}\right )}}
\]
[In]
integrate(sec(d*x+c)*(a+a*sec(d*x+c))^(5/2)*(A+C*sec(d*x+c)^2),x, algorithm="fricas")
[Out]
2/315*((903*A + 584*C)*a^2*cos(d*x + c)^4 + 2*(147*A + 146*C)*a^2*cos(d*x + c)^3 + 3*(21*A + 73*C)*a^2*cos(d*x
+ c)^2 + 130*C*a^2*cos(d*x + c) + 35*C*a^2)*sqrt((a*cos(d*x + c) + a)/cos(d*x + c))*sin(d*x + c)/(d*cos(d*x +
c)^5 + d*cos(d*x + c)^4)
Sympy [F]
\[
\int \sec (c+d x) (a+a \sec (c+d x))^{5/2} \left (A+C \sec ^2(c+d x)\right ) \, dx=\int \left (a \left (\sec {\left (c + d x \right )} + 1\right )\right )^{\frac {5}{2}} \left (A + C \sec ^{2}{\left (c + d x \right )}\right ) \sec {\left (c + d x \right )}\, dx
\]
[In]
integrate(sec(d*x+c)*(a+a*sec(d*x+c))**(5/2)*(A+C*sec(d*x+c)**2),x)
[Out]
Integral((a*(sec(c + d*x) + 1))**(5/2)*(A + C*sec(c + d*x)**2)*sec(c + d*x), x)
Maxima [F]
\[
\int \sec (c+d x) (a+a \sec (c+d x))^{5/2} \left (A+C \sec ^2(c+d x)\right ) \, dx=\int { {\left (C \sec \left (d x + c\right )^{2} + A\right )} {\left (a \sec \left (d x + c\right ) + a\right )}^{\frac {5}{2}} \sec \left (d x + c\right ) \,d x }
\]
[In]
integrate(sec(d*x+c)*(a+a*sec(d*x+c))^(5/2)*(A+C*sec(d*x+c)^2),x, algorithm="maxima")
[Out]
2/315*(315*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)^(1/4)*(5*(A*a^2*d*cos(2*d*x + 2*
c)^4 + A*a^2*d*sin(2*d*x + 2*c)^4 + 4*A*a^2*d*cos(2*d*x + 2*c)^3 + 6*A*a^2*d*cos(2*d*x + 2*c)^2 + 4*A*a^2*d*co
s(2*d*x + 2*c) + A*a^2*d + 2*(A*a^2*d*cos(2*d*x + 2*c)^2 + 2*A*a^2*d*cos(2*d*x + 2*c) + A*a^2*d)*sin(2*d*x + 2
*c)^2)*integrate((((cos(8*d*x + 8*c)*cos(2*d*x + 2*c) + 3*cos(6*d*x + 6*c)*cos(2*d*x + 2*c) + 3*cos(4*d*x + 4*
c)*cos(2*d*x + 2*c) + cos(2*d*x + 2*c)^2 + sin(8*d*x + 8*c)*sin(2*d*x + 2*c) + 3*sin(6*d*x + 6*c)*sin(2*d*x +
2*c) + 3*sin(4*d*x + 4*c)*sin(2*d*x + 2*c) + sin(2*d*x + 2*c)^2)*cos(9/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x +
2*c))) + (cos(2*d*x + 2*c)*sin(8*d*x + 8*c) + 3*cos(2*d*x + 2*c)*sin(6*d*x + 6*c) + 3*cos(2*d*x + 2*c)*sin(4*
d*x + 4*c) - cos(8*d*x + 8*c)*sin(2*d*x + 2*c) - 3*cos(6*d*x + 6*c)*sin(2*d*x + 2*c) - 3*cos(4*d*x + 4*c)*sin(
2*d*x + 2*c))*sin(9/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))))*cos(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*
x + 2*c) + 1)) - ((cos(2*d*x + 2*c)*sin(8*d*x + 8*c) + 3*cos(2*d*x + 2*c)*sin(6*d*x + 6*c) + 3*cos(2*d*x + 2*c
)*sin(4*d*x + 4*c) - cos(8*d*x + 8*c)*sin(2*d*x + 2*c) - 3*cos(6*d*x + 6*c)*sin(2*d*x + 2*c) - 3*cos(4*d*x + 4
*c)*sin(2*d*x + 2*c))*cos(9/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))) - (cos(8*d*x + 8*c)*cos(2*d*x + 2*c
) + 3*cos(6*d*x + 6*c)*cos(2*d*x + 2*c) + 3*cos(4*d*x + 4*c)*cos(2*d*x + 2*c) + cos(2*d*x + 2*c)^2 + sin(8*d*x
+ 8*c)*sin(2*d*x + 2*c) + 3*sin(6*d*x + 6*c)*sin(2*d*x + 2*c) + 3*sin(4*d*x + 4*c)*sin(2*d*x + 2*c) + sin(2*d
*x + 2*c)^2)*sin(9/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))))*sin(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x
+ 2*c) + 1)))/(((cos(2*d*x + 2*c)^4 + sin(2*d*x + 2*c)^4 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2
*d*x + 2*c) + 1)*cos(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos
(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c)^2 + 2*
cos(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(8*d*x + 8*c)^2 + 9
*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2
+ sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c)^2 + (2*cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c
) + 1)*sin(2*d*x + 2*c)^2 + 2*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^
2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2
+ 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(8*d*x + 8*c) + 6*(c
os(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2
*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(6*d*x + 6*c) + 6*(cos(2*d*x +
2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(4*d*x + 4*c) + co
s(2*d*x + 2*c)^2 + 2*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1
)*sin(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (
cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(8*d*x + 8*c) + 6*(sin(2*d*x + 2*c)^3 + 3*(c
os(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*co
s(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(6*d*x + 6*c) + 6*(sin(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + 2*cos(2
*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(4*d*x + 4*c))*cos(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))^
2 + (cos(2*d*x + 2*c)^4 + sin(2*d*x + 2*c)^4 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) +
1)*cos(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c)
^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c)^2 + 2*cos(2*d*x + 2
*c)^3 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(8*d*x + 8*c)^2 + 9*(cos(2*d*x +
2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x
+ 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c)^2 + (2*cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*
d*x + 2*c)^2 + 2*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x
+ 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*
x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(8*d*x + 8*c) + 6*(cos(2*d*x + 2*
c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) +
1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(6*d*x + 6*c) + 6*(cos(2*d*x + 2*c)^3 + cos
(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(4*d*x + 4*c) + cos(2*d*x + 2*c
)^2 + 2*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x +
6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2
*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(8*d*x + 8*c) + 6*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*
c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c
) + 1)*sin(2*d*x + 2*c))*sin(6*d*x + 6*c) + 6*(sin(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) +
1)*sin(2*d*x + 2*c))*sin(4*d*x + 4*c))*sin(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))^2)*(cos(2*d*x
+ 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)^(1/4)), x) + 20*((A + C)*a^2*d*cos(2*d*x + 2*c)^4 + (
A + C)*a^2*d*sin(2*d*x + 2*c)^4 + 4*(A + C)*a^2*d*cos(2*d*x + 2*c)^3 + 6*(A + C)*a^2*d*cos(2*d*x + 2*c)^2 + 4*
(A + C)*a^2*d*cos(2*d*x + 2*c) + (A + C)*a^2*d + 2*((A + C)*a^2*d*cos(2*d*x + 2*c)^2 + 2*(A + C)*a^2*d*cos(2*d
*x + 2*c) + (A + C)*a^2*d)*sin(2*d*x + 2*c)^2)*integrate((((cos(8*d*x + 8*c)*cos(2*d*x + 2*c) + 3*cos(6*d*x +
6*c)*cos(2*d*x + 2*c) + 3*cos(4*d*x + 4*c)*cos(2*d*x + 2*c) + cos(2*d*x + 2*c)^2 + sin(8*d*x + 8*c)*sin(2*d*x
+ 2*c) + 3*sin(6*d*x + 6*c)*sin(2*d*x + 2*c) + 3*sin(4*d*x + 4*c)*sin(2*d*x + 2*c) + sin(2*d*x + 2*c)^2)*cos(7
/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))) + (cos(2*d*x + 2*c)*sin(8*d*x + 8*c) + 3*cos(2*d*x + 2*c)*sin(
6*d*x + 6*c) + 3*cos(2*d*x + 2*c)*sin(4*d*x + 4*c) - cos(8*d*x + 8*c)*sin(2*d*x + 2*c) - 3*cos(6*d*x + 6*c)*si
n(2*d*x + 2*c) - 3*cos(4*d*x + 4*c)*sin(2*d*x + 2*c))*sin(7/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))))*co
s(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)) - ((cos(2*d*x + 2*c)*sin(8*d*x + 8*c) + 3*cos(2*d*x + 2
*c)*sin(6*d*x + 6*c) + 3*cos(2*d*x + 2*c)*sin(4*d*x + 4*c) - cos(8*d*x + 8*c)*sin(2*d*x + 2*c) - 3*cos(6*d*x +
6*c)*sin(2*d*x + 2*c) - 3*cos(4*d*x + 4*c)*sin(2*d*x + 2*c))*cos(7/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*
c))) - (cos(8*d*x + 8*c)*cos(2*d*x + 2*c) + 3*cos(6*d*x + 6*c)*cos(2*d*x + 2*c) + 3*cos(4*d*x + 4*c)*cos(2*d*x
+ 2*c) + cos(2*d*x + 2*c)^2 + sin(8*d*x + 8*c)*sin(2*d*x + 2*c) + 3*sin(6*d*x + 6*c)*sin(2*d*x + 2*c) + 3*sin
(4*d*x + 4*c)*sin(2*d*x + 2*c) + sin(2*d*x + 2*c)^2)*sin(7/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))))*sin
(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)))/(((cos(2*d*x + 2*c)^4 + sin(2*d*x + 2*c)^4 + (cos(2*d*x
+ 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*
x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2
*d*x + 2*c) + 1)*cos(4*d*x + 4*c)^2 + 2*cos(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(
2*d*x + 2*c) + 1)*sin(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*si
n(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c)^2 + (
2*cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c)^2 + 2*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*s
in(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c) + 3*
(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + c
os(2*d*x + 2*c))*cos(8*d*x + 8*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x
+ 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x +
2*c))*cos(6*d*x + 6*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)^2 +
cos(2*d*x + 2*c))*cos(4*d*x + 4*c) + cos(2*d*x + 2*c)^2 + 2*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin
(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*co
s(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(8*d
*x + 8*c) + 6*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4
*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(6*d*x + 6*c) + 6*(sin(2*d*x
+ 2*c)^3 + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(4*d*x + 4*c))*cos(5/2*arctan2(s
in(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))^2 + (cos(2*d*x + 2*c)^4 + sin(2*d*x + 2*c)^4 + (cos(2*d*x + 2*c)^2 + s
in(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 +
2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) +
1)*cos(4*d*x + 4*c)^2 + 2*cos(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)
+ 1)*sin(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c
)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c)^2 + (2*cos(2*d*x +
2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c)^2 + 2*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*
c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c) + 3*(cos(2*d*x +
2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*
c))*cos(8*d*x + 8*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + s
in(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(6*
d*x + 6*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x +
2*c))*cos(4*d*x + 4*c) + cos(2*d*x + 2*c)^2 + 2*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)
^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c
) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(8*d*x + 8*c) + 6
*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) +
(cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(6*d*x + 6*c) + 6*(sin(2*d*x + 2*c)^3 + (c
os(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(4*d*x + 4*c))*sin(5/2*arctan2(sin(2*d*x + 2*
c), cos(2*d*x + 2*c) + 1))^2)*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)^(1/4)), x) +
2*((13*A + 20*C)*a^2*d*cos(2*d*x + 2*c)^4 + (13*A + 20*C)*a^2*d*sin(2*d*x + 2*c)^4 + 4*(13*A + 20*C)*a^2*d*cos
(2*d*x + 2*c)^3 + 6*(13*A + 20*C)*a^2*d*cos(2*d*x + 2*c)^2 + 4*(13*A + 20*C)*a^2*d*cos(2*d*x + 2*c) + (13*A +
20*C)*a^2*d + 2*((13*A + 20*C)*a^2*d*cos(2*d*x + 2*c)^2 + 2*(13*A + 20*C)*a^2*d*cos(2*d*x + 2*c) + (13*A + 20*
C)*a^2*d)*sin(2*d*x + 2*c)^2)*integrate((((cos(8*d*x + 8*c)*cos(2*d*x + 2*c) + 3*cos(6*d*x + 6*c)*cos(2*d*x +
2*c) + 3*cos(4*d*x + 4*c)*cos(2*d*x + 2*c) + cos(2*d*x + 2*c)^2 + sin(8*d*x + 8*c)*sin(2*d*x + 2*c) + 3*sin(6*
d*x + 6*c)*sin(2*d*x + 2*c) + 3*sin(4*d*x + 4*c)*sin(2*d*x + 2*c) + sin(2*d*x + 2*c)^2)*cos(5/2*arctan2(sin(2*
d*x + 2*c), cos(2*d*x + 2*c))) + (cos(2*d*x + 2*c)*sin(8*d*x + 8*c) + 3*cos(2*d*x + 2*c)*sin(6*d*x + 6*c) + 3*
cos(2*d*x + 2*c)*sin(4*d*x + 4*c) - cos(8*d*x + 8*c)*sin(2*d*x + 2*c) - 3*cos(6*d*x + 6*c)*sin(2*d*x + 2*c) -
3*cos(4*d*x + 4*c)*sin(2*d*x + 2*c))*sin(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))))*cos(5/2*arctan2(sin
(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)) - ((cos(2*d*x + 2*c)*sin(8*d*x + 8*c) + 3*cos(2*d*x + 2*c)*sin(6*d*x + 6
*c) + 3*cos(2*d*x + 2*c)*sin(4*d*x + 4*c) - cos(8*d*x + 8*c)*sin(2*d*x + 2*c) - 3*cos(6*d*x + 6*c)*sin(2*d*x +
2*c) - 3*cos(4*d*x + 4*c)*sin(2*d*x + 2*c))*cos(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))) - (cos(8*d*x
+ 8*c)*cos(2*d*x + 2*c) + 3*cos(6*d*x + 6*c)*cos(2*d*x + 2*c) + 3*cos(4*d*x + 4*c)*cos(2*d*x + 2*c) + cos(2*d
*x + 2*c)^2 + sin(8*d*x + 8*c)*sin(2*d*x + 2*c) + 3*sin(6*d*x + 6*c)*sin(2*d*x + 2*c) + 3*sin(4*d*x + 4*c)*sin
(2*d*x + 2*c) + sin(2*d*x + 2*c)^2)*sin(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))))*sin(5/2*arctan2(sin(
2*d*x + 2*c), cos(2*d*x + 2*c) + 1)))/(((cos(2*d*x + 2*c)^4 + sin(2*d*x + 2*c)^4 + (cos(2*d*x + 2*c)^2 + sin(2
*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*co
s(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*
cos(4*d*x + 4*c)^2 + 2*cos(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)
*sin(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c)^2
+ 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c)^2 + (2*cos(2*d*x + 2*c
)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c)^2 + 2*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2
+ 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)
^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*
cos(8*d*x + 8*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2
*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(6*d*x
+ 6*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c)
)*cos(4*d*x + 4*c) + cos(2*d*x + 2*c)^2 + 2*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 +
2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) +
1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(8*d*x + 8*c) + 6*(si
n(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (co
s(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(6*d*x + 6*c) + 6*(sin(2*d*x + 2*c)^3 + (cos(2
*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(4*d*x + 4*c))*cos(5/2*arctan2(sin(2*d*x + 2*c),
cos(2*d*x + 2*c) + 1))^2 + (cos(2*d*x + 2*c)^4 + sin(2*d*x + 2*c)^4 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2
+ 2*cos(2*d*x + 2*c) + 1)*cos(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c
) + 1)*cos(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4
*c)^2 + 2*cos(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(8*d*x +
8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c)^2 + 9*(cos(2*d*
x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c)^2 + (2*cos(2*d*x + 2*c)^2 + 2*cos(2
*d*x + 2*c) + 1)*sin(2*d*x + 2*c)^2 + 2*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d
*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*
x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(8*d*x + 8
*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2
+ 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(6*d*x + 6*c) + 6*(c
os(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(4*d*x +
4*c) + cos(2*d*x + 2*c)^2 + 2*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x
+ 2*c) + 1)*sin(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x
+ 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(8*d*x + 8*c) + 6*(sin(2*d*x + 2*c
)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c
)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(6*d*x + 6*c) + 6*(sin(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2
+ 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(4*d*x + 4*c))*sin(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2
*c) + 1))^2)*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)^(1/4)), x) + 4*((3*A + C)*a^2*
d*cos(2*d*x + 2*c)^4 + (3*A + C)*a^2*d*sin(2*d*x + 2*c)^4 + 4*(3*A + C)*a^2*d*cos(2*d*x + 2*c)^3 + 6*(3*A + C)
*a^2*d*cos(2*d*x + 2*c)^2 + 4*(3*A + C)*a^2*d*cos(2*d*x + 2*c) + (3*A + C)*a^2*d + 2*((3*A + C)*a^2*d*cos(2*d*
x + 2*c)^2 + 2*(3*A + C)*a^2*d*cos(2*d*x + 2*c) + (3*A + C)*a^2*d)*sin(2*d*x + 2*c)^2)*integrate((((cos(8*d*x
+ 8*c)*cos(2*d*x + 2*c) + 3*cos(6*d*x + 6*c)*cos(2*d*x + 2*c) + 3*cos(4*d*x + 4*c)*cos(2*d*x + 2*c) + cos(2*d*
x + 2*c)^2 + sin(8*d*x + 8*c)*sin(2*d*x + 2*c) + 3*sin(6*d*x + 6*c)*sin(2*d*x + 2*c) + 3*sin(4*d*x + 4*c)*sin(
2*d*x + 2*c) + sin(2*d*x + 2*c)^2)*cos(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))) + (cos(2*d*x + 2*c)*si
n(8*d*x + 8*c) + 3*cos(2*d*x + 2*c)*sin(6*d*x + 6*c) + 3*cos(2*d*x + 2*c)*sin(4*d*x + 4*c) - cos(8*d*x + 8*c)*
sin(2*d*x + 2*c) - 3*cos(6*d*x + 6*c)*sin(2*d*x + 2*c) - 3*cos(4*d*x + 4*c)*sin(2*d*x + 2*c))*sin(3/2*arctan2(
sin(2*d*x + 2*c), cos(2*d*x + 2*c))))*cos(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)) - ((cos(2*d*x +
2*c)*sin(8*d*x + 8*c) + 3*cos(2*d*x + 2*c)*sin(6*d*x + 6*c) + 3*cos(2*d*x + 2*c)*sin(4*d*x + 4*c) - cos(8*d*x
+ 8*c)*sin(2*d*x + 2*c) - 3*cos(6*d*x + 6*c)*sin(2*d*x + 2*c) - 3*cos(4*d*x + 4*c)*sin(2*d*x + 2*c))*cos(3/2*
arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))) - (cos(8*d*x + 8*c)*cos(2*d*x + 2*c) + 3*cos(6*d*x + 6*c)*cos(2*d
*x + 2*c) + 3*cos(4*d*x + 4*c)*cos(2*d*x + 2*c) + cos(2*d*x + 2*c)^2 + sin(8*d*x + 8*c)*sin(2*d*x + 2*c) + 3*s
in(6*d*x + 6*c)*sin(2*d*x + 2*c) + 3*sin(4*d*x + 4*c)*sin(2*d*x + 2*c) + sin(2*d*x + 2*c)^2)*sin(3/2*arctan2(s
in(2*d*x + 2*c), cos(2*d*x + 2*c))))*sin(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)))/(((cos(2*d*x +
2*c)^4 + sin(2*d*x + 2*c)^4 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(8*d*x + 8
*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c)^2 + 9*(cos(2*d*x
+ 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c)^2 + 2*cos(2*d*x + 2*c)^3 + (cos(2*d*
x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d
*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(
2*d*x + 2*c) + 1)*sin(4*d*x + 4*c)^2 + (2*cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c)^2 + 2*
(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos
(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos
(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(8*d*x + 8*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x
+ 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4
*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(6*d*x + 6*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin
(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(4*d*x + 4*c) + cos(2*d*x + 2*c)^2 + 2*(sin(2*d*
x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c) + 3*(cos(2*
d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d
*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(8*d*x + 8*c) + 6*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x
+ 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x
+ 2*c))*sin(6*d*x + 6*c) + 6*(sin(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2
*c))*sin(4*d*x + 4*c))*cos(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))^2 + (cos(2*d*x + 2*c)^4 + sin(
2*d*x + 2*c)^4 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(8*d*x + 8*c)^2 + 9*(co
s(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + s
in(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c)^2 + 2*cos(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 +
sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 +
2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)
+ 1)*sin(4*d*x + 4*c)^2 + (2*cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c)^2 + 2*(cos(2*d*x +
2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)
+ 1)*cos(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c)
+ 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(8*d*x + 8*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*
d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2
*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(6*d*x + 6*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)
^2 + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(4*d*x + 4*c) + cos(2*d*x + 2*c)^2 + 2*(sin(2*d*x + 2*c)^3 +
3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2
+ sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1
)*sin(2*d*x + 2*c))*sin(8*d*x + 8*c) + 6*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*
cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(6
*d*x + 6*c) + 6*(sin(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(4*d*
x + 4*c))*sin(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))^2)*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2
+ 2*cos(2*d*x + 2*c) + 1)^(1/4)), x) + (A*a^2*d*cos(2*d*x + 2*c)^4 + A*a^2*d*sin(2*d*x + 2*c)^4 + 4*A*a^2*d*c
os(2*d*x + 2*c)^3 + 6*A*a^2*d*cos(2*d*x + 2*c)^2 + 4*A*a^2*d*cos(2*d*x + 2*c) + A*a^2*d + 2*(A*a^2*d*cos(2*d*x
+ 2*c)^2 + 2*A*a^2*d*cos(2*d*x + 2*c) + A*a^2*d)*sin(2*d*x + 2*c)^2)*integrate((((cos(8*d*x + 8*c)*cos(2*d*x
+ 2*c) + 3*cos(6*d*x + 6*c)*cos(2*d*x + 2*c) + 3*cos(4*d*x + 4*c)*cos(2*d*x + 2*c) + cos(2*d*x + 2*c)^2 + sin(
8*d*x + 8*c)*sin(2*d*x + 2*c) + 3*sin(6*d*x + 6*c)*sin(2*d*x + 2*c) + 3*sin(4*d*x + 4*c)*sin(2*d*x + 2*c) + si
n(2*d*x + 2*c)^2)*cos(1/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))) + (cos(2*d*x + 2*c)*sin(8*d*x + 8*c) +
3*cos(2*d*x + 2*c)*sin(6*d*x + 6*c) + 3*cos(2*d*x + 2*c)*sin(4*d*x + 4*c) - cos(8*d*x + 8*c)*sin(2*d*x + 2*c)
- 3*cos(6*d*x + 6*c)*sin(2*d*x + 2*c) - 3*cos(4*d*x + 4*c)*sin(2*d*x + 2*c))*sin(1/2*arctan2(sin(2*d*x + 2*c),
cos(2*d*x + 2*c))))*cos(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)) - ((cos(2*d*x + 2*c)*sin(8*d*x +
8*c) + 3*cos(2*d*x + 2*c)*sin(6*d*x + 6*c) + 3*cos(2*d*x + 2*c)*sin(4*d*x + 4*c) - cos(8*d*x + 8*c)*sin(2*d*x
+ 2*c) - 3*cos(6*d*x + 6*c)*sin(2*d*x + 2*c) - 3*cos(4*d*x + 4*c)*sin(2*d*x + 2*c))*cos(1/2*arctan2(sin(2*d*x
+ 2*c), cos(2*d*x + 2*c))) - (cos(8*d*x + 8*c)*cos(2*d*x + 2*c) + 3*cos(6*d*x + 6*c)*cos(2*d*x + 2*c) + 3*cos
(4*d*x + 4*c)*cos(2*d*x + 2*c) + cos(2*d*x + 2*c)^2 + sin(8*d*x + 8*c)*sin(2*d*x + 2*c) + 3*sin(6*d*x + 6*c)*s
in(2*d*x + 2*c) + 3*sin(4*d*x + 4*c)*sin(2*d*x + 2*c) + sin(2*d*x + 2*c)^2)*sin(1/2*arctan2(sin(2*d*x + 2*c),
cos(2*d*x + 2*c))))*sin(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)))/(((cos(2*d*x + 2*c)^4 + sin(2*d*
x + 2*c)^4 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(8*d*x + 8*c)^2 + 9*(cos(2*
d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2
*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c)^2 + 2*cos(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + sin(
2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*c
os(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)
*sin(4*d*x + 4*c)^2 + (2*cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c)^2 + 2*(cos(2*d*x + 2*c)
^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1
)*cos(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2
*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(8*d*x + 8*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x
+ 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x
+ 2*c)^2 + cos(2*d*x + 2*c))*cos(6*d*x + 6*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 +
2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(4*d*x + 4*c) + cos(2*d*x + 2*c)^2 + 2*(sin(2*d*x + 2*c)^3 + 3*(c
os(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + si
n(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*si
n(2*d*x + 2*c))*sin(8*d*x + 8*c) + 6*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(
2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(6*d*x
+ 6*c) + 6*(sin(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(4*d*x +
4*c))*cos(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))^2 + (cos(2*d*x + 2*c)^4 + sin(2*d*x + 2*c)^4 +
(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2
+ sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2
+ 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c)^2 + 2*cos(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^
2 + 2*cos(2*d*x + 2*c) + 1)*sin(8*d*x + 8*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*
c) + 1)*sin(6*d*x + 6*c)^2 + 9*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x +
4*c)^2 + (2*cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c)^2 + 2*(cos(2*d*x + 2*c)^3 + cos(2*d*
x + 2*c)*sin(2*d*x + 2*c)^2 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x +
6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x +
2*c)^2 + cos(2*d*x + 2*c))*cos(8*d*x + 8*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 3*
(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + c
os(2*d*x + 2*c))*cos(6*d*x + 6*c) + 6*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 2*cos(2*d*x
+ 2*c)^2 + cos(2*d*x + 2*c))*cos(4*d*x + 4*c) + cos(2*d*x + 2*c)^2 + 2*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*
c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c) + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c
)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c
))*sin(8*d*x + 8*c) + 6*(sin(2*d*x + 2*c)^3 + 3*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)
+ 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(6*d*x + 6*c) + 6*(
sin(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(4*d*x + 4*c))*sin(5/2
*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))^2)*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x +
2*c) + 1)^(1/4)), x))*sqrt(a) - (3*(105*A*a^2*sin(8*d*x + 8*c) + 140*(5*A + C)*a^2*sin(6*d*x + 6*c) + 126*(11*
A + 8*C)*a^2*sin(4*d*x + 4*c) + 12*(91*A + 73*C)*a^2*sin(2*d*x + 2*c))*cos(9/2*arctan2(sin(2*d*x + 2*c), cos(2
*d*x + 2*c) + 1)) - (315*A*a^2*cos(8*d*x + 8*c) + 420*(5*A + C)*a^2*cos(6*d*x + 6*c) + 378*(11*A + 8*C)*a^2*co
s(4*d*x + 4*c) + 36*(91*A + 73*C)*a^2*cos(2*d*x + 2*c) + (903*A + 584*C)*a^2)*sin(9/2*arctan2(sin(2*d*x + 2*c)
, cos(2*d*x + 2*c) + 1)))*sqrt(a))/((d*cos(2*d*x + 2*c)^4 + d*sin(2*d*x + 2*c)^4 + 4*d*cos(2*d*x + 2*c)^3 + 6*
d*cos(2*d*x + 2*c)^2 + 2*(d*cos(2*d*x + 2*c)^2 + 2*d*cos(2*d*x + 2*c) + d)*sin(2*d*x + 2*c)^2 + 4*d*cos(2*d*x
+ 2*c) + d)*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)^(1/4))
Giac [F]
\[
\int \sec (c+d x) (a+a \sec (c+d x))^{5/2} \left (A+C \sec ^2(c+d x)\right ) \, dx=\int { {\left (C \sec \left (d x + c\right )^{2} + A\right )} {\left (a \sec \left (d x + c\right ) + a\right )}^{\frac {5}{2}} \sec \left (d x + c\right ) \,d x }
\]
[In]
integrate(sec(d*x+c)*(a+a*sec(d*x+c))^(5/2)*(A+C*sec(d*x+c)^2),x, algorithm="giac")
[Out]
sage0*x
Mupad [B] (verification not implemented)
Time = 23.62 (sec) , antiderivative size = 766, normalized size of antiderivative = 4.53
\[
\int \sec (c+d x) (a+a \sec (c+d x))^{5/2} \left (A+C \sec ^2(c+d x)\right ) \, dx=\text {Too large to display}
\]
[In]
int(((A + C/cos(c + d*x)^2)*(a + a/cos(c + d*x))^(5/2))/cos(c + d*x),x)
[Out]
((a + a/(exp(- c*1i - d*x*1i)/2 + exp(c*1i + d*x*1i)/2))^(1/2)*((A*a^2*2i)/d - (a^2*exp(c*1i + d*x*1i)*(903*A
+ 584*C)*2i)/(315*d)))/(exp(c*1i + d*x*1i) + 1) - ((a + a/(exp(- c*1i - d*x*1i)/2 + exp(c*1i + d*x*1i)/2))^(1/
2)*(exp(c*1i + d*x*1i)*((A*a^2*2i)/d - (a^2*(A + 2*C)*4i)/d + (a^2*(21*A - 32*C)*2i)/(105*d)) - (A*a^2*2i)/(5*
d) + (a^2*(5*A + 2*C)*4i)/(5*d) - (a^2*(5*A + 32*C)*2i)/(5*d)))/((exp(c*1i + d*x*1i) + 1)*(exp(c*2i + d*x*2i)
+ 1)^2) - ((a + a/(exp(- c*1i - d*x*1i)/2 + exp(c*1i + d*x*1i)/2))^(1/2)*(exp(c*1i + d*x*1i)*((A*a^2*4i)/(3*d)
- (a^2*(3*A + C)*8i)/(9*d) + (a^2*(13*A + 20*C)*4i)/(9*d) - (a^2*(A + C)*40i)/(9*d)) - (A*a^2*4i)/(3*d) + (a^
2*(3*A + C)*8i)/(9*d) - (a^2*(13*A + 20*C)*4i)/(9*d) + (a^2*(A + C)*40i)/(9*d)))/((exp(c*1i + d*x*1i) + 1)*(ex
p(c*2i + d*x*2i) + 1)^4) + ((a + a/(exp(- c*1i - d*x*1i)/2 + exp(c*1i + d*x*1i)/2))^(1/2)*(exp(c*1i + d*x*1i)*
((A*a^2*10i)/(3*d) - (a^2*(189*A + 292*C)*2i)/(315*d)) - (A*a^2*2i)/(3*d) + (a^2*(9*A + 4*C)*2i)/(3*d)))/((exp
(c*1i + d*x*1i) + 1)*(exp(c*2i + d*x*2i) + 1)) + ((a + a/(exp(- c*1i - d*x*1i)/2 + exp(c*1i + d*x*1i)/2))^(1/2
)*(exp(c*1i + d*x*1i)*((A*a^2*10i)/(7*d) - (a^2*(3*A + 4*C)*10i)/(7*d) - (a^2*(9*A - 16*C)*2i)/(63*d) + (a^2*(
11*A + 20*C)*2i)/(7*d)) - (A*a^2*2i)/(7*d) + (a^2*(11*A + 4*C)*2i)/(7*d) - (a^2*(5*A + 12*C)*6i)/(7*d) + (a^2*
(5*A - 16*C)*2i)/(7*d)))/((exp(c*1i + d*x*1i) + 1)*(exp(c*2i + d*x*2i) + 1)^3)